The 15-20 fim emission in the reflection nebula NGC 2023 
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ABSTRACT 

We present 15-20 yum spectral maps towards the reflection nebula NGC 2023 ob- 
tained with the Infrared Spectrograph in short-wavelength, high-resolution mode on 
board the Spitzer Space Telescope. These spectra reveal emission from PAHs, Cdo, 
and H 2 superposed on a dust continuum. These emission components exhibit distinct 
spatial distributions: with increasing distance from the illuminating star, we observe 
the PAH emission followed by the dust continuum emission and the H 2 emission. The 
C 60 emission is located closest to the illuminating star in the south while in the north, 
it seems to be associated with the H/H 2 transition. Emission from PAHs and PAH- 
related species produce features at 15.8, 16.4, 17.4, and 17.8 yum and the 15-18 /mi 
plateau. These different PAH features show distinct spatial distributions. The 15.8 fim 
band and 15-18 //m plateau correlate with the 1 1 .2 yum PAH band and with each other, 
and are attributed to large, neutral PAHs. Conversely, the 16.4 /urn feature correlates 
with the 12.7 jxm PAH band, suggesting that both arise from species that are favored 
by the same conditions that favor PAH cations. The PAH contribution to the 17.4 fim 
band is displaced towards the illuminating star with respect to the 11.2 and 12.7 //m 
emission and is assigned to doubly ionized PAHs and/or a subset of cationic PAHs. 
The spatial distribution of the 17.8 yum band suggests it arises from both neutral and 
cationic PAHs. In contrast to their intensities, the profiles of the PAH bands and the 
15-18 //m plateau do not vary spatially. Consequently, we conclude that the carrier of 
the 15-18 /mi plateau is distinct from that of the PAH bands. 

Subject headings: Astrochemistry - Infrared : ISM - ISM : molecules - ISM : molecular data - 
ISM : line and bands - Line : identification - techniques : spectroscopy 
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1. Introduction 



Strong emission bands at 3.3, 6.2, 7.7 and 1 1.3 /mi - the so-called unidentified infrared (UIR) 
bands - dominate the mid-IR spectra of almost all objects, including reflection nebulae, planetary 
nebulae, the interstellar medium and HII regions. These bands are generally attributed to IR 
fluorescence of a family of Polycyclic Aromatic Hydrocarbon molecules (PAHs) pumped by the 
UV radiation field. A key result in the observational studies of PAHs is that their mid-IR bands 



show clear variations in peak positions, shapes an d (relative) intensities, not on 



but a 



2006 



so spatially within extended sources (e.g. 



Honvetal. 



Smith et al. 



2007b; 



Galliano et al. 



2001 



Peeters et al. 



y between sources 



2002; 



Brandl et al 



20081) . The variability in intensity ratio of the main 



PAH bands has been interpre t ed in terms of the c harge state and excitation levels of the PAHs 



(e.g. 



Allamandola et al 



1999; 



Galliano et al 



2008) wh ile the variation i n the profiles of the mai n 



PAH bands (classified as class A, B and C profiles; 



Peeters et al. 



2002 



2004) is attributed to the environment (CSM versus ISM; 



2002 



Peeters et al 



Boersma et al. 



Boersma et al 



2008). 



van Diedenhoven et al 



20021 



2008[) and the degree of UV processing (ISloan et al 



Van Kerckhoven 



2007 



Tielens 



2008; 



In addition to the well-studied main PAH bands, a plethora of weaker PAH bands are 
present, for example in the 15-20 /zm region. PAH emission in this region was first observed 
with the Infrared Space Observatory (ISO): new features were reported at 15.8, 16.4, 17.4 and 



17.9 urn with the 16.4 and then 17.4 /mi PAT 



1996 



Moutou et al. 



2000; 



Van Kerckhoven et al. 



Sturm et al. 



200ol : 



bands being the most 



Van Kerckhoven et al. 



promi nent (IBeintema et al. 



2000). In addition, 



(|2000 ) presented evidence for a variable, broad plateau from 15 to 20 um 



which seems to be present solely in HII regions (|Van Kerckhoven et al 



2000; 



from 15 to 2 


um 


Peeters et al. 


2004. 



2006). Based upon observations with the Spitzer Space Te 



are reported at 16.6, 17.0 and 17.2 jxm (ISmith et al. 



2004 



escope 



2007b; 



additional broad co mponents 



Sellgren et al. 



2007). But it 



is really the high sensitivity of Spitzer that revealed the omnipresence of these bands within 
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2004; 



Sellgren et al. 



2007 



Smith et al 



the Milky Way and in other galaxies (e.g. IWerner et al.ll2004bl: ISmith et al.ll2004l : iBrandl et al. 



2007b; 



Tappe et al 



2006), allowing - for the first time 



systematic investigations of these weaker bands. 



The 15-20 nm emission is attributed to C-C-C bending vibrations (e.g. l_A llamandola et al 



1989; 



Moutou et al. 



Boersma et al 



2010 



996; 



Van Kerckhoven et al. 



Ricca et al 



2000; 



Peeters et al. 



2004 



Mattioda et al. 



2009; 



2010Q . As a consequence, the bands in the 15-20 //m region are 



expected to reveal more about the overall molecular structure of the carriers than do the major 
bands below 15 /mi. An intriguing result of the Spitzer observations is the remarkable similarity 
(to first order) of the 15-20 /im PAH emission spectra from regions spanning a large range of 
physical conditions. This suggests that the astronomical PAH population may be dominated by a 
handful of stable, molecular structures (cf. survi val of the fittest) rath er than comprised of a large 



number of PAHs with widely varying structures (jBoersma et al 



20101 ). Hence, a systematic study 



of a large sample of spectra showing these weaker PAH bands will provide a complementary view 
on the characteristics of the emitting PAH population. 

This paper reports such a study. Here we analyze Spitzer- IRS spectral maps of NGC 2023 
in the 15-20 yum region. In Section |2j we describe the reflection nebula NGC 2023, while the 
observations and data reduction are discussed in Section [3] The data analysis is presented in 
Section |U and discussed in Section [5l We end with a short summary in Section [6] 



2. The reflection nebula NGC 2023 



NGC 2023, one of the best studied reflection nebuk 


le, is also among 


the 


brightest in the 


IR (e 




KnaDD et al. 1975 


; Harvey et a 


. 1980: 


Sellgren 


1984; 


Witt et al. 


1984 




Sellgren et al. 


1985; 


Gatlev et al. 1987: A 


Vitt & Ma 


in 


1989; I 


owe e 


tal 


. 1991 


Jaffe et al 


. 1994 


. Sellgren et al. 


1996; 


Steiman-Cameron et al. 1997; 


Field et al 


1998; 


Martini et all 1999: 


McCartnev etal. 1999: 
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Fig. 1.— The IRAC [8.0] image of NGC 2023 with the SH FOV shown (white) for both the north 
and south posi tions studied he re. The star HD37903 is indicated by a black circle, sources A, C 



and D are from 



Sellaren ( 



19831) . and the white circles indicate 2MASS point sources located inside 
the SH apertures. S refers to the south ridge, SSE to the south- southeastern ridge, N to the north 
ridge, and NW to the northwestern ridge. Maps shown in this paper use the orientation denoted by 
the vectors (X, Y). 



Mookeriea et alJl2 000 
2004bl : lsellgren et all 



; Takami et alJl2000h IWvrowski et alJl200ft Burgh et alJl2002t IWerner et al 
2007L boinh . Illuminated by the B1.5 V star HD 37903, it is at a distance 
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of 350 pc (IMookerjea et al. 



2009; 



ShefFer et al 



201 1|) . Images reveal a network o f narrow 



filam ents and hot spots (see e.g. Fig.Q]); the width of the H 2 filaments is ~0.25-0.5"(iRouaiLStaL 



1997|) . The amount of morphological detail in these images is particularly striking. While the 



overall morphology is similar at all wavelengths, there are significant variations in location and 
contrast between filaments and background in the different tracers. The limb-brightened shells 
so prominent in the various maps of this reflection nebula provide a unique opportunity to study 
the layered spatial distribution of Photo-Dissociation Regions (PDRs). Detailed analysis of the 



( Jaffe et al. 


1990; 


— ^ » 

Burton et al. 


1990; 



S teiman-Cameron et al 



1997 



Wvrowski et al 



Tielens & Hollenbach 



1997 



2000) 



1985b imply that a FUV 



in a few selected locations using PDR models (e.g. 
radiation field of 10 4 times the average interstellar value is incident on a clumpy molecular cloud, 



with the bulk of the gas having densities varying from 10 3 to 10 5 cm 3 , depending on location. 



3. The data 



3.1. Observations 



r 



The observations were taken with the Infrared Spectrograph (IRS, Houck et al. 



2004) 



2004a|) and were part of the Open Time 



on board the Spitzer Space Telescope (IWerner et al. 
Observations PIDs 20097 and 5051 1. 

We obtained spectral maps for two positions in the reflection nebula NGC 2023 (see Fig.Q]): 
towards the dense shell 78" south of the exciting star HD 37903 correspondi ng to the H2 emis sion 



1998). 



peak, and towards a region to the north (+33", +105") of the exciting star (Bur ton et al. 
The north position is characterized by a much lower density of ~ 10 4 cm" 3 tha n the density at 



the south position, which is > 10 5 cm 3 ([Burton et al. 



1998 



Sheffer et al. 



2011). Likewise, the 
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Table 1 : Observation log 





north position 


south position 


pi p\ 
r ID 


^fK1 1 
JUjI 1 


zuuy / 


AUK 


2633 /v24 


i4U33y2U 


cycles x ramp time 


2x30s 


1x30s 


pointings | 


7 


12 


step size | 


5" 


5.65" 


pointings _L 


15 


12 


step size _L 


2.3" 


4.7" 


background" 


5:42:1.00, -2:6:54.5 





a a, 6 (J2000); units of a are hours, minutes, and seconds, and 
units of 5 are degrees, arc minutes, and arc seconds. 



north position has a slightly lower UV radiatio n field between 6 and 13.6 eV, i.e. ~ 5 00 Gcu, as 



compared to ~ 10 4 G at the southern H 2 peak (|Burton et al.ll 1998c ISheffer et al. 



2011). 



The spectral maps are made with the short-high (SH) mode. The SH mode covers a 
wavelength range from 10 to 20 /mi at a resolution of ~ 600 and has a pixel size of 2.3". We 
obtained off-source background observations in SH for the north position. Table[T]gives a detailed 
overview of the observations. We also obtained spectral maps made with the short-low (SL) mode; 
the combined SL & SH data set will be presented in Peeters et al. (2012, in prep.). 



is the integrated 6 to 13.6 eV radiation flux in units of the Habing field = 1.6x10 3 
ergs/cm 2 /s. 
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3.2. Reduction 



The raw data were processed with the SI 8.7 pipeline versi on by the Spitzer S cience Center. 



2007a) and irsclean 



The resulting bed-products are further processed using cubism (|Smith et al. 
available from the SSC website^ 

The SH off-source background observation exhibits a rising dust continuum and an 1 1.2 /mi 
PAH band. Since we do not have off-source background observations for the south position, we 
investigated the importance of the background flux and, in particular, of the background PAH 
emission with respect to the on-source (PAH) flux. The 1 1.2 PAH flux in the SH background for 
the northern position is 8.5 10~ 8 W/m 2 /sr. This amounts to 3-11% of the on-source 11.2 PAH 
flux across the north map and to <1-10% of the on-source 1 1.2 PAH flux across the south map 



excluding source D from 



Sellgrenl (|1983|) . Source D is located on the edge of the south map 



and has a background contribution to the PAH flux of up to 45%. These results are consistent 
with the IRAC [8 /urn] image which is dominated by PAH emission and indicates a background 
flux of ~ 6% the minimum on-source flux. Hence, aside from source D, we conclude that the 
background has only a small contribution to the on-source PAH flux and therefore we did not 
apply a background subtraction. For the remaining of this paper, we therefor excluded source D 
from the PAH analysis. 

Since no background subtraction was applied, we used the tool irsclean to correct for the 
presence of rogue pixels. These rogue pixels are particularly problematic for the high-resolution 
module SH. In addition, we applied cubism's automatic bad pixel generation with (Ttrim = 7 and 
Minbad-fraction = 0.50 and 0.75 for the global bad pixels and record bad pixels respectively. 
Remaining bad pixels were subsequently removed manually. 



Spectra were extracted from the spectral maps by moving, in one pixel steps, a spectral 



2 http : //s sc . spitzer.caltech.edu 
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aperture of 2x2 pixels in both directions of the maps. This results in overlapping extraction 
apertures. No jumps in flux level were seen between the different orders within the SH module 
(SHI 1 to SH20) and hence no scaling was necessary. 



3.3. The spectra 



A typical 15-20 /mi spectrum of NGC2023 reveals a rising dust co ntinuum, a 



H 2 emission 



Sellgren et al. 



ine ( 



7 /mi), an emission band at 19 /mi attributed to C(,o (ICami et al 



sharp 



20101 : 



2010), a feature at 17.4 /mi which is a blend of a PAH band with another C^o 
band and a plethora of PAH emission bands at 15.8, 16.4, 17.4 and, 17.8 /mi on top of a broad 
underlying emission component (Fig. [2]). 



1983). The YSO source D is located 



NGC2023 contains a cluster of young stars (ISellgrenl . 
just outside the FOV of the south map (see Fig.Q])- Nevertheless, the spectra of regions close to 
source D show PAH emission features as well as characteristics typically for YSOs, i.e. a strong 
dust continuum (see the continuum map at 19.29 /urn in Fig. [3]) and a (strong) CO2 ice feature 
near 15 /mi. Furthermore, the background contribution to the 11.2 PAH flux is significant for a 
large fraction of these (ice-)spectra. We therefore excluded these spectra in the PAH and the 14.9 
/dm continuum analysis. The 15-20 //m spectrum of the YSO source C has the same spectral 
characteristics as the spectra across NGC 2023 but has an enhanced surface brightness (see pixels 
~ (21, 6) in the continuum maps at 14.9 and 19.29 /urn in Fig. [3]). However, we did not exclude 
these spectra in the analysis. 



3.4. Continuum and fluxes 



We determined the fluxes of the PAH features in the 15-20 //m region by subtracting a local 
spline continuum shown as a green line in the top frame in Fig. [21 This continuum is determined 
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by using anchor points at 14.91, 15.20, 15.50, 16.14, 16.69, 16.89, 17.16/17.19 (for the south 
and north maps, respectively), 17.58, 18.14/18.17 (for the north and south maps, respectively), 
18.50, 19.29 and 19.36 /urn. The fluxes are then estimated by fitting a Gaussian profile to the 
individual emission bands. A second global continuum (red line top frame) is determined by 
using anchor points at 14.91, 18.50 and 19.36 /mi. In this way, a plateau underneath the individual 
bands is defined by the difference of both continua and is further referred to as the 15-18 /im 
plateau. This plateau is disti nct from the 15-20 /im plateau observed towards HII regions 



(Van Kerckhoven et al 



2000j) which exhibits emission over the entire 15-20 /im range. This 



decomposition method is further referred to as the spline decomposition. 

The applied decomposition of the 15-20 /urn region is not unique. For example, 



Sellgren et al. 



(|2007j) fit this region by a combination of a continuum fit (similar to our global continuum) and 
eight Gaussians (further referred to as the G aussian decompositi on). Similarly, a combination 



of Drude profiles is regularly used (PAHFIT, 



Smith et al 



2007b|) . For comparison, we applied 



both methods to our spectra (Fig. [2]). To provide sufficiently wavelength coverage for PAHFIT, 
we obtained the SL spectrum for the same extraction aperture as the SH spectrum. However, 
we did not correct for the different PSF at these wavelengths. The PAHFIT decomposition in 
the 15-20 jjm region results in components repr esenting the dust continuum emission (which 



is a combination of modified blackbodies, see 



Smith et al. 



2007bl) . the H 2 emission and the 



PAH dust features including a broad co mponent centered at 17.0 /im. We also applied a similar 



decomposition as 



Sellgren et al 



(|2007L their fig. 8 and table 1) to our spectra by using their 
Gaussian parameters as a starting value. This method results in two possible decompositions of 
the 15-20 /im region; one dominated by a broad component peaking around ~ 17.2 /im (Fig. 
[21 third panel) or another one in which this 17.2 /mi component is much narrower resulting in 
broader and stronger components near 16.6 and 17.4 /mi (Fig. [21 fourth panel). Hence, in order 
to study the PAH bands in a consistent way, one needs to restrict the decomposition to a broad or 
narrow 17.2 /im component. But choosing one over the other will clearly influence the derived 
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fluxes of the Gaussian components. 

The Gaussian decomposition with the broad 17.2 /mi component results in very similar 
fluxes for the features as those obtained with the spline decomposition discussed above, except for 
the 15.8 and 16.4 //m bands (Fig. [2l bottom panel). In contrast, different PAH band intensities are 
obtained with the Gaussian decomposition with the narrow 17.2 /urn component and the PAHFIT 
decomposition. In particular, the broad component centered around 17 /jm depends significantly 
on the decomposition method. Here we will adopt the spline decomposition (Fig. [21 top panel). 

4. Data analysis 

Here we investigate i) the relationship between individual PAH emission bands and ii) the 
relationship between the different emission components - PAHs, H 2 , C6o, and dust continuum 
- present in the 15-20 /mi region based upon the spatial distribution and the intensity ratios of 
the various emission components. We further investigate the PAH band profiles in the 15-20 /mi 
region. The implications of these results are discussed in the next section. 

4.1. Relationship between emission components 

Maps: Figs. |3]and|4]show the spatial distribution of the various emission components in the 
15-20 /mi region. In these figures, in addition to the color-coded maps for the different features, 
we show contours for the 11.2 and 12.7 /mi PAH bands for reference. These were derived from 
the SH data and are discussed in Peeters et al. (2012, in prep.). The north map is characterized by 
lower flux levels and therefore larger scatter is present in the map tracing the 15.8 and 17.8 //m 
features and, to a lesser extent, the 18.9 /mi band. Hence, when discussing the spatial distribution 
of the 15.8 and 17.8 /mi bands, we restrict ourselves to the south map. 
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Feature correlations: To investigate possible variations within the PAH spectrum itself, 
we need to exclude the influence of PAH abundance and column density. This is achieved by 
normalizing the fluxes to another band in this region. Here we start with the 11.2 /mi PAH band 



which is attributed to solo CH out-of -plane bending vibration of large, neutral PAHs (IHony et al. 



2001 



Bauschlicher et al 



2008, 



20091) . To look for possible correlations with the 1 1.2 /mi band, 



we also apply a normalization to the 12.7 /mi PAH band which is attributed to duo and trio CI 



out-of-plane bending vibration of large PAHs (|Hony et al 



2001 



Bauschlicher et al. 



2008, 



2009) 



Both the 1 1.2 and 1 2.7 /mi PAH flux es are obtained from the SH data using a spline continuum 



(see e.g. Fig. 1 of IHony et al 



2001) and are discussed in Peeters et al. (2012, in prep.). A 
correlation is often found between the normalized intensity ratios of PAH bands in this region as 
these bands arise from a family of PAH molecules which all have the same basic structural unit. 
However, variations in charge states, excitation level, edge structures, etc., exist within the PAH 
family which can tighten as well as loosen correlations between the normalized intensity ratios 
of specific PAH bands. For example, the well known tight correlation betwe en the 6.2 and 7.7 



um PAH bands is a ttributed to the common charge state of these bands (e.g. 



Galliano et al. 



Hony et al. 



2001 



2008J). Figs. [5] and [6] show some of the intensity correlations found in NGC 2023. 



The correlation coefficients and the fit parameters are given in Tabled While maps smooth out the 
detailed, fine grained variations, allowing one to see the big picture, these correlations highlight 
subtle differences. 

Results: The following trends are derived from the south spectral maps shown in Fig. [3j 
The 11.2 /mi feature, the 15.8 /mi feature and the 15-18 /mi plateau show very similar spatial 
morphology with distinct peaks at the S and SSE positions in the IRAC map. In contrast, the 
distribution of the 12.7, 16.4, 17.4, and 17.8 //m features are displaced towards the star. The 12.7 
and 16.4 /mi bands have very similar spatial behavior; both peak at the SSE position and show a 
broad, diffuse plateau NW of the line connecting the S and SSE ridges. The 17.4 //m feature also 
presents a strong, diffuse plateau NW of that line, but does not show a peak at the SSE position 
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as does the 12.7 and 16.4 /zm bands. In addition, the 17.4 /im feature has a separate peak in the 
upper right of the field, coinciding with the peak in the 18.9 /im band. The 17.8 /im band peaks 
primarily on the SSE position, overlapping with the peak of both the 1 1.2 and 12.7 /im bands but 
also shows some enhanced emission slightly offset from the 1 1.2 /zm peak near the S ridge, though 
not as much as the 12.7 and 16.4 /zm bands. The 18.9 /mi band is clearly in a class by itself, 
peaking closer to the illuminating star than any of the other features considered here. Interestingly, 
a second peak is evident in the 18.9 /im spatial distribution at source C. Note also that the 18.9 
/im band does show (weak) emission tracing the S and SSE ridges. The spatial behavior of the 
14.9 and 19.29 /zm continuum intensities align closely with the H 2 emission distribution, all lying 
along the SE portion of the 11.2 /im contours, while clearly avoiding regions with strong 12.7 /zm 
emission. The H 2 map shows, however, additional emission towards the west. 

The similarity in distribution of the 11.2 /im and 15-18 /zm plateau is also evident in the 
north map (Fig. H]): both maps peak on the NW ridge and extend to include the N ridge. There is 
too much scatter in the weak 15.8 /zm PAH band to warrant comparison. In addition, the 16.4 /zm 
PAH band has a spatial distribution very similar to that of the 12.7 /zm band. However, unlike 
the situation at the SSE ridge where the 11.2, 12.7 and 16.4 /zm bands all peak in the same place, 
the 12.7 and 16.4 /zm bands are clearly displaced towards the south with respect to the 11.2 /zm 
emission in both the peak emission and the N ridge extension. The 17.4 and 18.9 /zm bands 
exhibit more scatter in the north map. Nevertheless, the 17.4 /zm band seems to peak towards the 
NW ridge and further west of this ridge. In contrast, the 18.9 /zm band seems to trace the N ridge 
and not the NW ridge while also showing secondary emission in the south. The spatial behavior 
of the 14.9 and 19.29 /zm continuum emission and the H 2 emissions are quite distinct from the 
PAH emission in the north map as they do not track the 1 1.2 /zm contours as closely and only 
trace the N ridge and not the NW ridge. 

This morphology suggests that the features can be viewed as falling into the following 
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groups: 1) the 11.2, 15.8 //m bands and the 15-18 /mi plateau; 2) the 12.7, and 16.4 /urn bands; 
3) the 17.4 /mi diffuse component; 4) the 17.8 //m band; 5) the 18.9 and (some of) the 17.4 /mi 
features; 6) the underlying continuum producing the emission measured here at 14.9 and 19.2 /mi; 
and 7) the H 2 emission. 

The correlation plots shown in Fig. [5] support grouping the 11.2, 15.8 //m bands and the 
15-18 /mi plateau into a single group. Indeed, a very tight correlation is found in our data set 
between the 11.2/12.7 PAH intensity ratio and the 15-18 /mi plateau/12.7 PAH intensity ratio. 
Weaker trends, although still very tight, are seen between the 11.2/12.7 and 15.8/12.7 PAH 
intensity ratios and the 15.8/12.7 and 15-18 //m plateau/12.7 PAH intensity ratios. In addition, the 
normalized 16.4 //m PAH intensity clearly correlates with that of the 12.7 /mi PAH band; their 
correlation coefficient is the highest found in our data, i.e. 0.889. 

As shown in Fig. [6j the 17.4 and 18.9 //m bands do not exhibit an overriding linear 
correlation with each other as do the bands in Fig. [5] Instead, the 18.9 /mi band does not correlate 
with the 17.4 /urn band for low 18.9/1 1.2 intensity ratios (< 0.02) and shows a linear correlation 
with the 17.4 //m band for strong band ratios. This is consistent with the assignment of the 18.9 
urn bandto Cfip and the 17.4 /mi feat ure to being a blend of a PAH emission band with a C 60 band 



(ICami et al. 



20101 : 



Sellgren et al. 



2010|) . Therefore, the 17.4 //m band was decomposed into these 
two components to estimate the spatial distribution of the PAH contribution to this band. This is 
referred to as the diffuse 17.4 /mi band. We calculated the average ratio of the 18.9/17.4 band 
intensity in the upper right corner of the south map (where the 18.9 /urn band peaks) and used this 
value to subtract the Ceo contribution to the 17.4 /mi band to obtain the "diffuse 17.4 //m band" 
for both the south and north map (see Fig. [3] and®. Note that this method assumes a constant 
ratio between the 17.4 and 18.9 //m C6o emission (of 0.6, similar to what is observed e.g. in Tel 



(Camietal 



20101) ) and only provides a first order approximation to the diffuse 17.4 /mi band 



distribution. Nevertheless, this lends support to our earlier suggestion that the PAH component in 
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the 17.4 fxm feature does not spatially coincide with the 11.2 or 12.7 /mi PAH emission. 

4.2. Band profiles 

We investigate the band profiles in the spectrum from the total FOV of the south map to 
minimize the noise and apply the same continuum determination as discussed in Sect. 13.41 Within 
the uncertainty, the 15.8, 17.4, 17.8, and 18.9 //m bands are well fitted with a Gaussian. Hence, 
no shift is found in the central wavelength of the 17.4 /im feature as the 17.4/18.9 intensity ratio 
varies. In contrast, the 16.4 /im PAH feature is slightly asymmetric to the red (Fig. [7]). To 
assess the influence of the continuum on the degree of asymmetry, three different continua were 
constructed and the resulting 16.4 /im bands compared. First, the continuum point near 16.69 
/im may be considered to be part of the red wing of the 16.4 /mi band. Hence, it is omitted 
and a spline is fitted to the remainding continuum points. The other two ways involve fitting a 
straight line to continuum points at 16.14 /mi and 16.69 or 16.89 /mi respectively. There is little 
difference between the spline fit and the straight line but the omission of the continuum point at 
16.69 //m clearly amplifies the red wing of the band and thus the asymmetry (Fig. [7]). Comparison 
of the 16.4 /mi band profile with the well-known asymmetric 6.2 and 11.2 /mi PAH bands shows 
that its asymmetry is less pronounced (Fig. [7]). 

In order to investigate possible variations in the profile of the 15-18 /mi plateau, we compared 
the 15-18 /mi plateau of each pixel with the average 15-18 /mi plateau (Fig. [8]). Overall, we 
conclude that the profile of the 15-18 /im plateau does not vary, despite a possible change in 
intensity of the PAH bands located on top of this plateau. Only in a handful of cases do we see 
slight deviation from the average. Hence, the 15-18 /mi plateau profile does not depend on the 
individual PAH bands as determined in this paper. 
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5. Discussion 



5.1. Correlation studies 



In a previous de tailed study of the PAH emission features in the 15-20 //m region, 



Boersma et al. 



(|201Qf) focused on the spectra from very different astronomical environments, not 



the spatial behavior in an e xtend ed object as considered here. Analyzing the spectra from some 



10 objects, 



Boersma et al. 



(2010) showed that the 16.4 fxm band correlated with the mid-IR 6.2 



and 7.6/7.8 fim bands and that the total 15-20 //m PAH emission correlated with the 11.2 //m 
band. No further connections were found between the 15.8, 16.4, 17.4, 17.8 and 18.9 /mi bands 
or between the spectral characteristics of the 15-20 //m emission and the mid-IR PAH modes. 
These authors then conclude that these bands must be carried by independent molecular species 
or classes of species. The results presented in the previous section both support a nd amplify t hese 



Boersma et al. 



conclu sions. The strong correlation of the 16.4 and the 6.2 //m band found by 
(|2010|) is confirmed by our results as both the 16.4 and the 12.7 /mi PAH bands strongly correlate 
with each other a nd it is well known that the 12.7 /mi correlate s well with the 6.2 /urn band 



(|Hony et al.l l2.Q01. Peeters et al. 2012, in prep.). In contrast with lBoersma et al. 



(120101) . we find 



a correlation between the 15.8 and 11.2 /mi bands. The 15.8 /mi band is very weak and only 
detected in a few sources in their sample. These sources follow the correlation found in this paper. 



5.2. The 15-18 /urn plateau 



As discussed in Sect. 13 .4[ different spectral decompositions have been applied in the 
literature. This is in particular true for the treatment of the broad underlying plateaus. They are 



either considered to b e part of the PAI 



or Drude profiles (e.g. 



Boulanger et al. 



bands themselves which are then fitted by Lorentzians 



1998 



Smith et al 



independent of these PAH bands (e.g. iHony et al 



2001; 



2007b), or t hey are considered to be 



Peeters et al. 



2002, this paper). The 
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spatial maps and the band correlations clearly indicate that there are significant differences in the 
spatial distribution between the 15-18 pm plateau and most PAH bands. These spatial differences 
suggest a different carrier for the plateau and the individual PAH bands (with the exclusion of 



the 15.8 pm band), a finding consistent with ear lier suggestions for t 



prominent mid-IR bands between 5 and 13 pm (|Bregman et al 



19891 : 



le plateaus underl ying the 



Roche et al 



I989L 



5.3. The 15-20 pm emission components 



The emission features between 15-20 /mi can be collected into seven groups as discussed in 
Section |4?T1 Here we discuss these in terms of the specific characteristics that set each group apart 
and make them particularly sensitive to the local physical conditions. 

Group 1: The 11.2, 15.8 /mi bands and the 15-18 pm plateau. Since the 11.2 um band is 



assoc i ated with the solo CH out-of-plane vibration of large, neutral PAHs (e 



1999 



Hony et al 



2001 



Galliano et al. 



2008 



Bauschlicher et al. 



2008, 



g- 



Allamandola et al. 



20091), the 15.8 /mi band 



and the 15-18 pm plateau must arise from species that are favored by the same conditions that 
favor this subset of the PAH population. In contrast with expectations, the 15-20 pm region 
shows little systematic dependence on P AH class or molecular structure as present in the NASA 



Ame s PAH IR spectroscopic database (|Bauschlicher et al. 



20101 : 



Boersma et al. 



2010! : 



Ricca et al. 



2010Q . suggesting that molecular structure is not the main origin for the observed correlations. 
Keep in mind that, to date, the database is biased to smaller PAHs and the astronomical emission 
at these longer wavelengths is dominated by larger PAHs. Thus, as knowledge of larger PAH 
spectroscopy increases, this suggestion should be revisited. In ad dition, ionization only has a 



moderate influence o n the PAH emission in the 15-20 pm region (|Van Kerckhoven et al 



Boersma et al. 



2000; 



20101) : it slightly influences the relative strength of the emission bands but not their 



peak position. Based on theoretical and experimental PAH spectra, one can therefore not assign 
a particular charge state to these individual bands, unlike for the mid-IR bands. Therefore, we 
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suggest that large, neutral PAHs are also responsible for the 15.8 /im band and 15-18 /im plateau. 
PAH clusters, VSGs (Very Small Grains) or small HAC particles may also contribute to the 15-18 
/im plateau. However, the spread of the correlations in Fig. [5] shows there is not a one-to-one 
correspondence between the carriers of these different emission features. These differences reflect 
minor changes in the relative populations of the different band carriers, perhaps arising from 
slightly different molecular structures, excitation levels, and so on. 

Group 2: The 12.7 and 16.4 /im bands. Similarly, since the 16.4 and 12.7 um b ands correlate 



with the 6.2 and 7.7 //m PAH cation bands (Ho ny et al J 12001c LBoersma et al 



2010, Peeters et 



al. 2012, in prep.), the 12.7 and 16.4 /im bands must arise from species that are co-spatial with 
the cationic portion of the emitting PAH population. Again, as with Group 1, the spread of the 
correlations in Fig. [5] shows there is not a one-to-one correspondence between these different 
carriers. However, in contrast with Group 1, we cannot propose molecular structure or charge as 
the main origin for these bands. Indeed, laboratory and theoretical PAH spectra indicate that the 
12.7 and 16.4 /im bands are sensitive to the molecular edge structure: the former being due to 



duo and trio CH groups (Ho ny et alJl2001c iBauschlic her et al. 



systematically present in PAHs with pendent r ings (IMoutou et alj 



2000 



Peeters et al. 



2004; 



Boersma et al. 



2008 



199 



2009) and th e latter being 

I 



Van Kerckhoven et al. 



2010 ) and large PAHs with pointy edges (|Ricca et al. 



20101) . Laboratory and theoretical PAH spectra al so indicate that both neutral and ionized 



can contribute to the 12.7 and 1 6.4 urn emission (IVan Kerckhoven et al. 



B au schlicher et al 



2008 



2009; 



Boersma et al 



2000; 



Hony et al 



AHs 



2001 



2010 ) and hence do not allow, by themselves, the 



assignment of these bands to a single charge state. Similarly, blind signal separation of the PAH 
emission across PDRs cannot assign a single charge state to the 12.7 jxm band (the 16.4 /im band 
is not analysed by this method). Indeed, this method revealed three spatially different components, 
referred to as PAH , PAH + an d VSG, which are attributed to respectively neutral PAHs, ioni zed 



PAHs, and PAH clusters (e.g. 



Rapacioli et al 



2005 



Berne et al. 



2007 



Rosenberg et al 



201 lh . 



The 12.7 fim band is present in both the PAH and PAH + components. However, the relative 
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contribution of the 12.7 //m band in the PAH° and PAH + components seems t o vary for differen t 



decomposition methods and sources (IRapacioli et al. 



2005 



Berne et al 



2007; 



Joblin et al 



2008) 



Hence, further investigation is warranted. 



Group 3: The 17.4 ixm diffuse component. This paper confirms the results by Sellgren et al. 
(120101) that the 17.4 pm emission band is composed of emiss ion due to C 60 and PAHs (i.e., the 



diffuse 17.4 /im component). In addition, 



Sellgren et al. 



(120 10|) report a similar spatial distribution 



between the diffuse 17.4 //m emission and the 16.4 /im PAH band in the reflection nebulae 
NGC7023. In contrast, in NGC 2023, the spatial distribution of the diffuse 17.4 pm emission is 
quite distinct from that of the 16.4 pm band. In particular, in the north map, the diffuse 17.4 pm 
emission is strong close to the peak 16.4 pm emission as well as west of the NW ridge. In the 
south map, it is located closer towards the illuminating star compared to the 1 1.2, 12.7 and 16.4 



2.7 and 16.4 um 



AH 



20101) 



fxm PAH emission and does not peak in the SSE ridge where the 1 1.2, 
bands peak. Perusal of the NASA Ames PAH IR spectroscopic database (|Bauschlicher et al. 
reveals that large compact PAHs systematically show emission near 17.4 pm and that PAH charge 
only has moderate influence on the relative intensities a nd no influence on the peak position of 



the bands in the 15-20 pm region (|Boersma et al. 



2010). We therefore conclude that the diffuse 



17.4 pm band is due to doubly ionized PAHs and/or a subset of PAH cations that favor more 
harsh conditions compared to those responsible for the 6.2 /im PAH bands, such as for example 
dehydrogenated cations. This possible assignment will further be explored by comparison to the 
PAH emission at shorter wavelengths in a forthcoming publication (Peeters et al. 2012, in prep.). 

Group 4: The 17.8 pm band. The 17.8 pm emission exhibits spatial characteristics between 
those of the 11.2 and 12.7 pm bands. With the former spatial distribution being attributed to 
neutral PAHs and the latter being co-spatial with that of the PAH cations, the 17.8 /im band likely 
originates from both neutrals and cations (with both having similar intrinsic strengths for this 
band). 
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Group 5: T he 18.9 and (some of) the 17.4 ixm fea ture. These two bands have recently been 



assigned to C 60 (iCami et al. 



2010 



Sellgren et al. 



2Q10Q . The C^q bands show a distinct spatial 



distribution. In the south, they reach their peak emission intensities in regions closer to the excitin g 



star than are the PAH bands, similar to what is observed towards NGC7023 (ISellgren et al. 



2010) 



In addition, secondary peaks in emission are found near source C, a YSO, and where the PAH 
emission peaks. In contrast, in the north, the peak in the Cgo distribution seems to be farther away 
from the illuminating source (no other UV sources are present near this position) and linked to 
the H/H 2 transition zone. This different behavior in the C6o emission may well provide important 
clues to the formation and evolution of Cgo in the ISM. This will be further explored when the full 
PAH spectra will be analyzed (Peeters et al. 2012, in prep.). 

Group 6: The continuum producing the emission at 14.9 and 19.2 pm. This emission 
component peaks further into the filament with respect to the PAH emission. This is consistent 
with the results of the blind signal separation dis cussed above: the VSG component, attri buted to 



PAH clusters, carries the mid-IR continuum (e.g. 



Rapacioli et al 



2005; 



Berne et al 



2007). As it is 



located deeper in the PDR compared to the PAH emission, these authors suggested the destruction 
of PAH clusters by UV photons and subsequent formation of PAHs. As mentioned above, PAH 
clusters may also contribute to the 15-18 pm plateau, which is slightly offset compared to the dust 
continuum emission. This suggest that the carriers of these emission components slightly differ. 



G roup 7: The H? emiss ion. 



refer to 



Fleming et al 



(120101) and 



or the analysis of t he H 2 emission present in this data set, we 



Sheffer et al 



(2011) 



6. Conclusions 



We studied the emission in the 15-20 pm region towards the reflection nebula NGC 2023 
by analysing two spectral maps obtained with Spitzer's Infrared Spectrograph (IRS), short- 
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wavelength high-resolution mode. We observed PAH emission bands at 15.8, 16.4, 17.4, and 17.8 
/mi, a broad PAH plateau between 15-18 /mi, C 60 emission at 17.4 and 18.9 /mi, and H 2 emission 
at 17.0 /mi superposed on a dust continuum. 



We found distinct spatial distributions for these emission components and found tight 
intensity correlations between some emission components. Based on these results, we collected 
the emission components in seven groups and discussed their specific characteristics. 

Group 1: the 11.2, 15.8 /mi PAH bands and the 15-18 /urn PAH plateau. We attributed this 
group to large, neutral PAHs. 

Group 2: the 12. 7 and 16.4 /mi PAH bands. Compared to Group 1 , these bands are displaced 
towards the illuminating star. We concluded that they must arise from species that are co-spatial 
with the cationic portion of the emitting PAH population. 

Group 3: the 17.4 /mi diffuse component. We estimated the PAH contribution to the 17.4 //m 
band and found that its spatial distribution is displaced towards the illuminating star compared 
to Groups 1 and 2. We assigned this PAH band to doubly ionized PAHs and/or a subset of the 
cationic PAH population, such as for example dehydrogenated PAHs. 

Group 4: the 17.8 //m PAH band. This band has a spatial distribution in between that of 
Groups 1 and 2. We suggested this band arises from both neutral and ionised PAHs. 

Group 5: the 18.9 and (some of) the 17.4 /mi features. The C fin emission shows a sp atial 



Sellgren et al 



(2010). 



distribution distinct from that of the PAH emission, consistent with 
Moreover, the spatial behavior of Cgo towards NGC2023 is different in the south and north. 
Specifically, the Cgo emission is located closer to the illuminating star compared to the PAH and 
H 2 emission in the south while in the north, it seems to be associated to the H/H 2 transition. 

Group 6: the underlying continuum producing the emission measured here at 14.9 and 
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19.2 /mi. The continuum emission is strongest deeper into the PDR compared to the PAH 
emission. This i s cons iste nt with the VSG co mponent attributed to PAH clusters by, for example, 



Rapacioli et al. 



(2005J) and 



Berne et al. 



(120071) 



Group 7: the H 2 emission. The H 2 emission is located furthest into the PDR. 



We also investigated the profiles of the PAH bands and 15-18 /im PAH plateau. All PAH 
bands have a symmetric band profile except for the 16.4 //m band, which is slightly asymmetric. 
We did not observe any changes in the band profiles spatially across the maps. Likewise, we found 
that the profile of the 15-18 jim PAH plateau does not vary. We concluded that the carrier of the 
underlying 15-18 /urn plateau is distinct from the individual PAH bands located on top of it. 

EP owes her genuine gratitude to Martin Houde for his help and support throughout this 
project. We would like to thank the referee Dr. K. Sellgren whose comments have helped to 
improve the paper. We very gratefully acknowledge sustained support from the NASA Spitzer 
Space Telescope General Observer Program. Studies of interstellar PAHs at Leiden Observatory 
are supported through advanced-ERC grant 246976 from the European Research Council. 
LJA gratefully acknowledges sustained support from NASA's Laboratory Astrophysics and 
Astrobiology Programs. 
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Fig. 2. — Typical 15-20 /urn spectra towards NGC 2023. The top panel exemplifies the spline de- 
composition applied in this paper: the local continuum (green), global continuum (red) and the 15- 
18 /um plateau (lined region). For comparison, the decomposition obtained with PAHFIT is shown 
in the second panel and the decomposition obtained by applying the Gaussian decomposition used 



by 



Sellgren et al. 



(2007) for their analysis of NGC7023 is shown in the middle panel. This Gaus- 
sian decomposition is also applied to a second spectrum towards NGC 2023 (fourth panel). The 
colour coding is as follows: continuum (red), fit (green), individual components (blue), the broad 
17 /um component (PAHFIT, purple) and the broad 17.2 yum component (Gaussian decomposition, 
magenta). The bottom panel shows the residual spectra obtained from the observed spectrum 1 by 
subtracting the local spline continuum, the 17 yum PAHFIT component and both the 16.6 and 17.2 
yum Gaussian components respectively. 




Fig. 3. — Spatial distribution of the emission features in the 15-20 /jm region in NGC 2023 for the south map. As a reference, 
the intensity profiles of the 11.2 and 12.7 fjm emission features (Peeters et al. 2012, in prep.) are shown as contours in black 
and white, respectively. The white arrow in the top right corners indicates the direction towards the central star and the E-N 
orientation is given in the bottom right panel (see Fig.Q]). S refers to the south ridge, SSE to the south-southeastern ridge, C to 
source C and D to source D. The axis labels refer to pixel numbers. The PAH fluxes and the 14.9 jim continuum flux in regions 
near source D excluded from the analysis are set to zero (near pixel (10,0); see Section [33T) . 
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Fig. 4. — Spatial distribution of the emission features in the 15-20 yum region in NGC 2023 for the north map. As a reference, 
the intensity profiles of the 11.2 and 12.7 fjm emission features (Peeters et al. 2012, in prep.) are shown as contours in black 
and white, respectively. The white arrow in the top right corners indicates the direction towards the central star and the E-N 
orientation is given in the bottom right panel (see Fig. [I])- N refers to the north ridge and NW to the northwestern ridge. The 
axis labels refer to pixel numbers. 
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Fig. 5. — Correlations in PAH intensity ratios across NGC 2023. The red squares are for the south 
map and the blue triangles for the north map. The fit to the data is shown as a solid black line and 
its parameters can be found in Table [2] The green line represents a fit to the data forced through 
(0,0). 
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Fig. 6. — The 17.4 versus the 18.9 /im band intensities across NGC2023. The red squares repre- 
sent the south map and the blue triangles the north map. The fit to the data for high band intensity 
ratios is shown as a black line. 

Table 2: Parameters of the fit to the observed correlations and their Pearson's correlation coeffi- 
cients. 



PAH intensity ratios 


correlation 


A 1 


B 1 


y vs. x 


coefficient 






15.8/12.7 vs. 11.2/12.7 


0.834 


-0.02 


0.02 


plat./12.7 2 vs. 11.2/12.7 


0.867 


0.21 


0.26 


plat./12.7 2 vs. 15.8/12.7 


0.837 


0.56 


8.79 


16.4/11.2 vs. 12.7/11.2 


0.889 


0.02 


0.16 



Parameters of the fit: y = A + Bx 
2 Plat, refers to the 15-18 //m plateau. 
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Fig. 7. — The asymmetric 16.4 /mi band (top panel). Two profiles are shown to assess the influence 
of the continuum determination on band shape: a spline fit to the continuum points discussed in 
Sect. 13.41 (black solid line) and a spline fit to the same continuum points excluding that at 16.69 
yum (grey solid line). The asymmetry of the black profile is assessed in the second top panel: the 
dotted red line shows the red wing mirrored along a vertical position at the peak position and the 
dashed blue line represents the mirrored blue wing. For reference, the asymmetric 6.2 and 11.2 fim 
band profiles are shown in the bottom and middle panel respectively. Note that the 6.2 band profile 
is obtained from the IRS-SL observations and the 1 1 .2 band profile from the IRS-SH observations. 
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Wavelength (fJ-m) 

Fig. 8. — The 15-18 /im plateau. Shown are spectra from different pixels (black), the average 
15-18 /im plateau scaled such that its flux matches the 15-18 /urn plateau flux of the black spectra 
(blue) and, the residuals (green). The second panel shows the extend of the variation seen in the 
15-18 jum plateau profile. Note that the strength of the features varies independently of the 15-18 
fim plateau and of each other. 
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